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We present preliminary results of the calculation of the positive-parity ground state charm baryon 
spectrum using of A?/ = 2 + 1 + 1 dynamical quarks. The calculation uses a relativistic heavy- 
quark action for the valence charm quark, clover- Wilson fermions for the valence light quarks 
and HISQ sea quarks. The spectrum is calculated with a lightest pion mass around 220 MeV 
and two lattice spacings {a ^ 0.12 fm and 0.09 fm) are used to extrapolate to continuum limit. 
Our preliminary results are consistent with the currently measured baryon spectrum, except for 
the isospin-averaged 7=1/2 'Ecc which is approximately 2 a above the SELEX observed value. 
We predict the yet-to-be-discovered double and triple-charm baryons E*^., £lcc, ^cc^ ^ccc to have 
masses 3665(42)(29) MeV, 3694(40)(45) MeV, 3739(35)(21) MeV and 4782(24)(28) MeV, re- 
spectively. 
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1. Introduction 

In recent years, interest in charmed baryon spectroscopy has resurfaced. This excitement 
has been partly triggered by the first observation of a double charm baryon (SELEX) candidate 
S+ (3520) [1] as well as its isospin partner S++(3460) [2]. The SELEX Collaboration has later 
confirmed their observation of S+ (3520) [3], but the BABAR [4], BELLE [5], and FOCUS [6] 
experiments have seen no evidence for either state of the isospin doublet (S+,E++). The SELEX 
evidence for this doublet implies unprecedented dynamics, the most surprising of which is the 
60 MeV mass difference between the two states. All other previously observed isospin partners 
have mass differences one order of magnitude smaller. A nice discussion of the theoretical status 
of the double charm spectrum can be found in [7], one remarkable feature is the fact that previous 
unquenched Lattice Quantum Chromodynamics (LQCD) calculations suggest the Ecc mass to be 
approximately 100-200 MeV higher than the one observed by the SELEX collaboration [8, 9]. 

There remain many undiscovered doubly and triply charmed baryons states. The recently up- 
graded Beijing Electron Positron Collider (BEPCII) detector, Beijing Spectrometer (BES-III), the 
LHC, and the GSI future project, the antiProton ANnihilation at DArmstadt (PANDA) experiment, 
will help further disentangle the heavy-baryon spectrum and resolve puzzles like the one mentioned 
above. LQCD calculations serve as direct first-principles theoretical input for these experiments. 
With this goal in mind, we present preliminary results of the first LQCD calculation of the charm 
baryon spectrum using Nf = 2 + 1 + 1 dynamical quarks. 

2. Lattice Formulation and Details 

In this work we use the gauge configurations generated by the MILC Collaboration [10], in 
which the Nf = 2+1 + 1 highly improved staggered quark (HISQ) action is used for the sea quarks. 
The implementation of the HISQ action has shown to further reduce lattice artifacts as compared 
to asqtad action[10]. The gauge configurations are HYP-Smeared to remove ultraviolet noise from 
the gauge field. The valence light-quark (up, down and strange) propagators are generated using 
the clover- Wilson action. 

In order to systematically remove the large ^((mcC?)") discretization artifacts (where nic is the 
charm-quai^k mass) we use the following implementation of the relativistic heavy-quark action for 
the charm sector [11-14]: 

Sq = L {m + YoDo -\dI + v (jiDi - ^of) - ^CBOijFij - ^CEao,Fo/)^, Q',, (2. 1) 



where Qx is the heavy-quark field at the site x, 7v are the Hermitian gamma matrices that satisfy the 
Euclidean Clifford algebra, a^v = ^[Y^iiTv], is the first-order lattice derivative, and F^v is the 
Yang-Mills field strength tensor. The parameters {mo, v,cb,ce} must be tuned to assure i^{{mca)") 
terms have been removed. For the coefficients cb and ce we use the tree-level tadpole-improved 
results [8, 15] cb = ^,c-e = with the tadpole factor, uq, defined as mq = I (LpTr(?7p))^/"*. 

The coefficients mo and V are simultaneously determined nonperturbatively by requiring the 
ratio ^ = '"'^"^J^^''^'* to be equal to its experimental value, 0.084, and {tj^ , to satisfy the cor- 
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masses cited above include the statistical and fitting window errors. Listed are the two bare masses of the 
valence charm-quarks used, the number of configurations, and the number of measurements performed for 
each ensemble. 

rect dispersion relation, E]^ = mjj + p^. The dispersion relation is matched using rjc and en- 
ergies for the at the six lowest momenta: ^(0,0,0), ^(1,0,0), ^(1,1,0), ^(1,1,1), ^(2,0,0), 
^(2,1,0), and their permutations. Two charm-quark masses are used for each ensemble. Fur- 
thermore, in order to account for further discretization systematic errors, we have performed the 
calculation at two lattice spacings, a ~ (.09 fm,0.12 fm); at each lattice spacing we use two dif- 
ferent light-quark masses corresponding to mji ~ 220 MeV,310 MeV. This is all done in order 
to be able to extrapolate to the physical point = 0.007, mj^/mjj = 0.084,(3 = fm). The 

lattice spacing for the ensembles is determined by extrapolating aniQ linearly in m\/m^ to its 
physical value. Our course (fine) lattice spacings are within 1.5 a (1 a) from the values obtained 
by the MILC collaboration [10] and within 1.7 a (1 a) from the determination by the HPQCD 
collaboration using the T 25 — 15 splitting [16]. 

Details of the ensembles, including our determination of the lattice spacing are listed in TA- 
BLE 1 . The systematic uncertainty in the pion mass presented in TABLE 1 due to the placement 
of the fitting window is calculated as the standard deviation of all fitting windows within three 
timeslices of the window picked for performing the fit. Although calculations are performed at a 
single volume, from the l^mji values listed in TABLE 1 it is clear that these calculations fall within 
the p-regime, ensuring that finite volume artifacts are exponentially suppressed. 

The interpolating operators used for the positive-parity ground state charmed baryons are the 
following [17]: 



J=l/2 


J=3/2 


A, = e^^'-P+Q'l {q[jT^q'2) 

i'L.y = e'^^'-P+iPl'^y^Q'^f, {q'jTiq'-) 
{KV = '-^PHPe^'Y^Q^- {q'JrJqf+qi^rJq:) 
i^,y = e'""P+{P'J^yjQ', {qi^jq^) 

{Z,,y = e^^'-P+{Pl'^)'iql (Gf n'G?) 
{accV = e'''"'P+{Py^yiq'l [Q^TiQ^) 


{Ky = 'iip^ipp'rQ'c WJrJqf+q'/rJq'::) 

{a*y = e''^"'P+{Pl'^)'iQ'l {q^jTjqf) 

{KcV = e'^^p^ipT^ql (ef n'G?) 

{Q.l,y = e'^"-P+{Pl'^yjq'l (Gf HG^") 



where quj^s respectively denote the up, down and strange quark annihilation operators, Qc de- 
notes the charm-quark operator, {k,l,m} are color indices, while depict polarization indices. 
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(r'*,P) are the antisymmetric and symmetric spin matrices (Cj^jCY), where C is the charge con- 
jugation matrix. In order to have the best possible overlap with the state of interest, we have used 



the spin projection operators {P^' )ij = 5ij - ^YiYj and {P^' )ij = 5ij - {P^' )ij = \yiyj, and the 
positive-parity projection operator = (1 + 74) /2. 
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Figure 1: Sample effective mass plots from the various ensembles. The fit includes two errorbars, the 
smallest of which is the statistical uncertainty, while the wider band includes the statistical and system- 
atic uncertainty added in quadrature. Additionally listed are P^r degree of freedom and the Q{d) = 
2dPY{d/2) -^X^ '^XoiXoY ^^^'^ factor which depends on the number of degrees of freedom d and is opti- 

mally near 1 . 

All two-point functions were calculated with gauge-invariant Gaussian-smeared sources and 
point sinks. Four sources are used in each gauge configuration.The masses were extracted from the 
single exponential fit to the correlation function at large Euclidean time. Double-exponential form 
is also used to crosscheck the ground-state mass. The uncertainties of the extracted hadron masses 
are evaluated using the jackknife method. Figure (1) displays examples of the effective mass plots 
for various coiTclation functions. The plots show the per degree of freedom, the Q(d) factor, 
along with the statistical and fitting window uncertainties for the fitting window chosen. 

3. Results and Conclusion 

For every hadron we calculated the ratio of its mass to the Q. mass, niH/niQ,. The mass of 
each hadron is calculated at two different values of the charm-quark mass, two pion masses, and 
two lattice spacings. The first of these allows us to interpolate the ratio, niH/ma, linearly in anic 
to the physical charm-quark mass defined by ^ = 0.084. We then simultaneously extrapolated 
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Figure 2: Shown is a survey of previous lattice calculations [8, 9, 18-21], along with the results of this 
paper labeled as "Briceno et al." The statistical uncertainty is show as a thick/inner error bar, while the 
statistical and systematic uncertainties added in quadrature are shown as a larger thin error bar. Our system- 
atic uncertainties include errors originating from the fitting window and scale setting. The experimentally 
determined masses are listed for comparison [22]. It is clear that our results are in agreement with the ex- 
perimental values (< a) with the clear exception of Z^c, where our result is 1.9c7's above the experimental 
value. 



the four values of the hadron masses. To perform the extrapolation we used the leading order 
polynomial dependance on mji and on the lattice spacing, a, mn/mci = ro + c\m\/m^-\-C2a^- The 
continuum limit hadron mass is recovered by multiplying the physical Q. mass, 1672.45(49) MeV. 
Using this procedure, we have verified that our calculations reproduce the low-lying charmed me- 
son spectrum. We check on the residual heavy-quark discretization effects using the mesonic 
spectrum after extrapolating to the physical point, and we obtain the following results for the 
charmonium hyperfine splitting, D and Dj masses: Ai^ = rtiji^i — nij^ = 111.0(5. 5)(6.6) MeV, 
mu = 1862.5(7.9)(6.7) MeV, mD,=1965(6)(ll) MeV. These are well in agreement with their cor- 
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responding experimentally observed values (116 MeV, 1869.6 MeV, 1968.5 MeV respectively). 
The hyperfine splitting demonstrates the necessity for continuum extrapolation when evaluating 
the charmed spectrum. The values of our hyperfine splitting for the course and fine ensembles 
are A%^^ = 76(12)(8) MeV, A?^"'' = 95.0(4.5)(4.9) MeV, the finest of which underestimates the 
physical value by approximately 18%. This underestimation is consistent with previous lattice 
calculations [23, 24]. 

In Fig. (2) we display our preliminary charmed baryon results, along with a survey of previous 
lattice calculation of the J=l/2 and J=3/2 ground state spectrum and the corresponding experimen- 
tal values. All previous calculations of the charmed baryon spectrum have been performed with 
Ught-quark masses corresponding to > 290 MeV, placing our calculations closest to the phys- 
ical point. Additionally, from Fig. (2) it is clear that our calculation represents the most precise 
determination of the continuum spectrum of charmed baryons. Our results are well in agreement 
with the experimental masses for all particles observed with the clear exception of E^c, where our 
result for the courses lattice is 1.9c7's above the experimental value. This ~ 120 MeV deviation 
from the experimentally determined mass has been observed in previous unquenched lattice calcu- 
lations [8, 9]. As mentioned above, the experimental value of this hadron has only been observed 
by SELEX collaboration with no confirmation by other groups. Theoretical prediction, our results 
included, suggests that this experimental finding by the SELEX collaboration needs further con- 
firmation. Lastly, we predict the yet-discovered double and triple-charm baryons masses S*^, D.cc, 
Q*,, Qccc to be 3665(42)(29) MeV, 3694(40)(45) MeV, 3739(35)(21) MeV and 4782(24)(28) MeV, 
respectively. 
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